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Lithium intercalated anatase used in Li-ion batteries has some special features: coexistence of Li-rich and Li- 
poor phases as well as two possible positions for Li ions in the oxygen tetrahedron. A theoretical description 
of the compound considering those peculiarities is presented. As shown by the performed symmetry analysis, 
the intercalation induced lattice deformation can be accompanied by the ordering of antiferroelectric type 
(internal piezoeffect). In the following step, a qualitative illustration of the phase separation in the lithiated 
anatase is given within the Landau expansion at the proper choice of coefficients. A microscopic model for 
description of the compound is also proposed which combines features of the Mitsui and Blume-Emery-Griffits 
models and utilizes the symmetry analysis results. Various ground state and temperature-dependent phase 
diagrams of the model are studied to find a set of model parameters corresponding to the lithiated anatase. 
A phase separation into the empty and half-filled phases in a wide temperature range has been found closely 
resembling the phase coexistence in the intercalated crystal. In the framework of the model, the two-position 
Li subsystem could have the ordering of ferro- or antiferroelectric types which, however, has not been yet 
observed by the experiment. 
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1. Titania in a nutshell 

1.1. Titanium dioxide: polymorphs, properties and applications 

There are eleven known polymorphs of titanium dioxide (titania) . The most common natural 
forms (rutile, anatase and brookite) are just different space arrangements of the TiC>6 group where 
a titanium cation is located in the centre of the slightly deformed octahedron shaped by six oxygen 
anions. Each polymorph has its own deviations from the ideal octahedron but it is always elongated 
along the certain axis. Thus, two (apical) oxygens are further away from the titanium than four 
others (equatorial). 

Being the most stable polymorphs, rutile and anatase are widely used and intensively studied. 
They are very similar in many details (e.g. the arrangement of atoms and average lengths of 
bonds) [ Ij. However, anatase is 10% less dense than rutile and an additional volume is condensed 
in the voids affecting the cell averaged properties such as compressibility and dielectric constant. 
Moreover, this minor difference becomes crucial at crystal intercalation. In both polymorphs, an 
elementary cell consists of two formula units but unlike the tetragonal rutile (P4 2 /mnra) the 
standard crystallographic cell of body-centered tetragonal anatase (74i / amd) is chosen to consist 
of two elementary cells. 

The titanium atoms, and hence, the octahedra, are arranged in such a way that each oxygen 
is at the same time an equatorial atom for one titanium, and an apical one for the other titanium 
atom in the same unit cell. Neighboring octahedra are sharing edges and corners with each other. 
Two and four edges of each octahedron are shared in rutile and anatase, respectively. The basic 
octahedra are distorted in such a way that each shared edge is shortened, the other edges being 
correspondingly elongated. The shortened oxygen-oxygen bonds are often named as the bridge 
bonds (in the sense that it bridges the interaction between Ti ions: metal-oxygen- metal) . In rutile 
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the bridge bond connects two equatorial oxygen atoms. Hence, the octahedra form the vertical 
linear chains. The octahedra belonging to adjacent chains are connected only through one corner: 
an oxygen atom which is both apical and equatorial for the two touching octahedra. In anatase, 
the octahedra are arranged in order to share a diagonal edge between an apical and an equatorial 
atom. Thus, octahedra form zig-zag chains orthogonal to the crystallographic axis. There are two 
sets of chains orthogonal to each other, that are connected through a common octahedron. 

The list of hi-tech applications of titanium dioxide is quite impressive. They are primarily 
related to its photoactivity. For example, TiC>2, particularly in the anatase form, is a photocatalyst 
under ultraviolet light. The strong oxidative potential of the positive holes oxidizes water to create 
hydroxyl radicals (the Honda- Fujishima effect [[2]). 

Superhydrophilicity phenomenon for glass coated with titanium dioxide is caused by ultraviolet 
light partially removing oxygen atoms from the surface of the titanium oxide. The areas where 
oxygen atoms were removed became hydrophilic, while the same size areas where no oxygen atoms 
were taken away turned out to be hydrophobic [[3|. The result is a TiOvcoated glass which is 
antifogging and self-cleaning. 

Dielectric properties of titanium dioxide distinguish it as semiconductor [ 0] (to stress the 
difference: anatase is semiconductor of the n-type while rutile is of p-type which is utilized in the 
gas sensor [;5]). Due to its high dielectric constant, it is commonly used as a dielectric in electronic 
devices, such as thin film capacitors [[6] and MOS devices [[7], as well as for the fabrication of 
anti-reflection coatings, interference filters [[8], as well as optical wave-guides [[9]. 

Although both rutile and anatase are potentially interesting for photo-catalysis and photo- 
electrochemical applications, experimental investigations have mostly focused on the more prospec- 
tive anatase polymorph. It has a wider optical-absorption gap and a smaller electron effective mass 
which presumably leads to a higher mobility for the charge carriers [ 10J and plays a key role in the 
injection process of novel dye-sensitized photochemical solar cells with high conversion efficiency [ 
lllj . In some cases, the two materials are used together in the same device, exploiting their peculiar 
properties for different purposes. For example, a typical low-cost photo- voltaic module is composed 
of a transparent conducting photo-electrode of dye-sensitized nanocrystalline anatase, a spacer of 
electrically insulating, light-reflecting particles of rutile, and a counter-electrode of graphite pow- 
der [[H|. Hence, anatase is used due to its efficient coupling with the dye, and rutile for its high 
dielectric constant. 

Furthermore, the open crystallographic structure of anatase facilitates the accommodation of 
substantial amounts of small ions (Li, H, etc.) within the lattice. Lithium insertion changes the 
optical properties of Ti02: it turns the white powder dark blue whilst in thin film form it changes 
from being transparent to partially reflecting (electrochromism) which is used in displays and sun- 
blinds (switchable mirrors) . Combining an electrochromic film and a photovoltaic film to form the 
two electrodes of an electrochemical cell one can achieve a photochromic structure [I13j. 

1.2. Intercalation of rutile and anatase: experiment and theory 

In the last decades the Li-ion batteries have run into operation as a result of their high energy 
capacity, re-chargeability and environmentally friendly properties. Anatase TiC>2 may act as an an- 
ode in such a battery [tH]. In practice, anatase is not the ideal candidate because of its relatively 
low potential versus other electrode materials. Better properties are demonstrated, e.g. by similar 
compounds with the spinel structure LiTi 2 04 [[H], Li4Ti 5 0i2 [HI] and their manganesian ana- 
logues LiMn 2 04 and LiMgo.1Nio.4Mn1.5O4 [HI]. Here anatase is considered as a well defined model 
material displaying many typical properties of transition metal oxide electrodes. The electrode and 
electrochromic properties of lithiated anatase are already well documented and partly exploited 
commercially. However, the impressive experimental breakthrough in the study of the microscopic 
processes resulting in these achievements is not accompanied by theoretical investigations. 

As mentioned earlier, anatase TiC-2 has a body-centered tetragonal structure indexed by the 
Iii/amd space group. Upon lithiation, anatase lattice undergoes an orthorhombic distortion that 
results in the Lio.sTiO^ phase (sometimes referred to as Li-titanate) indexed by the space group 
Imma, where the fourth order axis is lost due to the distortion in the ab plane [[18]. The change 
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in symmetry is accompanied by a decrease of the unit cell along the c-axis and by an increase 
along the 6-axis, resulting in a 4% increase of the unit cell volume. Lithium was found to reside in 
the interstitial voids within the oxygen octahedra [US]. The structural change can be explained as 
occupation of Ti-Ti bonding atomic orbitals by the electron that enters the TiC>2 lattice with each 
Li-ion to maintain charge neutrality. 

In comparison with the number of experimental studies of intercalated titanium dioxide the 
list of theoretical works on the subject looks very short containing primarily ab initio approaches. 
Some experimentally established properties are fairly well explained but some predictions are not 
confirmed by experiment. 

Already in the pioneer work [ [19] a higher possibility of lithium intercalation in the anatase 
structure than in rutile was predicted as well as the absorption energies obtained were calculated 
and Li-induced local one-electron energy levels were found in the gap between the upper valence 
band and the conduction band and could be attributed to Ti 3+ states. 

The calculations of the relative lithium insertion potentials were performed for the rutile, 
anatase, brookite, ramsdellite, colombite, spinel, and orthorhombic polymorphs of titanium diox- 
ide from the first principles periodic Hartree-Fock approach [[20] also indicating that lithium was 
completely ionized in LiTiC>2 and that the charge transfer is predominantly to the oxygen sublat- 
tice. A similar study of the average voltage to intercalate lithium in various metal oxides (among 
them TiC>2) and dichalcogenides was performed utilizing the ab initio pseudopotential method [ 
121] . It was also found that Li was fully ionized in the intercalated compounds but with its charge 
distributed among the anion and the metal. 

A series of works [ [22J [23l [231 ESI, [26], [23. [28] should be mentioned, where the lithium intercalation 
in both rutile and anatase was ab initio modelled taking into account thermodynamic and kinetic 
effects. The important role of strong local deformations of the lattice and elastic screening of 
interlithium interactions was established, the absence of insertion into rutile at room temperature 
was explained in terms of inaccessibility of the low-energy configurations due to highly anisotropic 
diffusion, a phase separation in anatase into a Li-rich phase and a Li-poor phase was described and 
the existence of a new distorted rock-salt phase for LiTiC>2 was predicted [[22l [28]. The calculated 
open circuit voltage profile reproduced and explained the characteristic features of experimental 
discharge curves for both polymorphs [[24]. An analysis of the site preference for Li intercalation 
in rutile and diffusion pathways of ions was performed. The expansion of the host structure on 
Li insertion was found to contribute to the enhanced diffusion of Li ions along the c direction 
while a large distortion of the rutile framework nearly suppressed Li diffusion in the ab planes; 
computed diffusion coefficients were found in excellent agreement with the measured values [ 26J. 
A new phase of LiTiC>2 is predicted which may be accessed through electrochemical lithiation of 
ramsdellite-structured TiC>2 at the lowest potential (remaining constant over a wide range of Li 
concentrations) reported for titanium dioxide based materials [ 29j. 

Ab initio study of the elastic properties of single and polycrystal TiC>2 and other IV-B group 
oxides in the cotunnite structure was performed [130]. 

Up till now only one non-"a& initio" description of the intercalation in titanium dioxide by means 
of the pseudospin-electron model [[31] (where the pseudospin formalism was used in describing the 
intercalant subsystem) was performed. The possibility of the phase transitions accompanied by an 
abrupt change of the concentration of intercalated ions and a significant increase of electrostatic 
capacity of the system was predicted. 

1.3. Special features of the lithium intercalated anatase: phase equilibrium and double posi- 
tions for lithium ions 

Upon lithium insertion, an increasing fraction of the material changes its crystallographic struc- 
ture from anatase TiC>2 to Li-rich lithium titanate Li .6TiO2 (sometimes a different stoichiometry 
is reported: Li .5TiO2 or Li .5sTiO2; as will be shown below it depends on the size of TiC>2 crys- 
tallites). Phase separation occurs on the Li-rich and the Li-poor (Lio.oiTiC^) phases [132], 

Such a two-phase equilibrium system in the electrodes provides a constant electrical poten- 
tial between their electrodes (so-called plateau in potential on the discharge curve) for a wide 
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range of the lithium concentration, because only the relative phase fractions vary on charging (or 
discharging) the lithium while their stoichiometries remain unchanged [[33]. The Li-rich lithium 
titanate phase progressively moves inside the anatase electrode as a front parallel to the interface 
and returns during lithium extraction exactly in the way it came in [ 34|. 

Both in the anatase and in the lithium titanate lattice, Li is found to be hopping over the 
available sites with activation energies of 0.2 and 0.09 eV, respectively. However, macroscopic 
intercalation data show activation energies of 0.5 eV because the diffusion through the phase 
boundary determines the activation energy of the overall diffusion and the overall diffusion rate 
itself [[32]. 

Recent NMR spectroscopy study [[35] of nanosized lithiated anatase revealed further important 
details of the phase behavior and morphology. The coexistence of the Li-poor and the Li- rich phases 
is possible only in the particles of the size exceeding 120 nm due to the surface strain (occurring 
between the phases) which becomes energetically unfavourable in small particles. For the system 
of 40 nm particles, phase stoichiometries are not stable indicating an enhanced solid solution 
behavior: lithium content increases to Li .iTiO2 in Li-poor and to Li .7TiO2 in Li-rich phases. 
Further decrease of the particle size makes it possible to find a fully occupied phase LiiTiC>2 [ 36j. 
It can coexist with the Li-rich one penetrating to the 3-4 nm depth below the particle surface (or 
transforming all the particles less than 7 nm in size). The poor Li ion conductivity can be due to 
the full occupation of the octahedral voids, whereas ion diffusion requires vacancies. Most likely, 
the short diffusion path in nanostructured materials diminishes this problem as well as elevation 
of temperature. 

As established by quasi-elastic neutron scattering [ 37J Li ions can occupy two distinct positions 
within the octahedral interstices along the c axis (but only one of them at a time) . In the Li-anatase 
those positions are symmetrical, separated by 1.61 A and equally occupied while in the Li-titanate 
they are shifted, separated by only 0.7 A and nonequivalent (riLii = 0.32 and nu2 — 0.19 at 10 K 
- the fitting of the site occupancy assuming a Boltzmann distribution indicates that the energy 
difference between the positions is 3.8 meV; positions 1 and 2 have an antiparallel orientation in 
the neighbouring octahedra due to the phase symmetry). A combination of quasi-elastic neutron 
scattering and force field molecular dynamics simulations shows that Li is hopping on a picosecond 
time scale between the two sites in the octahedral interstices [l38"ll39|. 

Lithium was also found to occupy multiple positions inside the distorted oxygen octahedron of 
LlrMgo.1Nio.4Mn1.5O4 spinel [[T7]. Quite possible that this feature is common to a wide family of 
crystals with a similar structure but has been found only recently due to a higher precision of the 
experiment. 



2. Symmetry analysis of the lithium intercalated anatase: a possibility of 
internal piezoeffect 



As mentioned before, the standard crystal- 
lographic cell of body-centered anatase is cho- 
sen to consist of two elementary cells for con- 
venience. Thus, the respective space group D 1 ^ 
(or IAi/amd origin choice 2, No. 141) has a 
double set of symmetry operations as compared 
with its point symmetry group D^h (^/mmm). 
However, a single elementary cell obeys all 
transformation rules. It consists of two formula 
units, i.e. two octahedral voids (formed by oxy- 
gen anions) where the intercalated lithium ion 
can reside in one of two available positions (fig- 
ure [1]). 

Occupation of each octahedron by the 
lithium ion could be easily described by the 
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Figure 1. A schematic illustration of the posi- 
tions available for the intercalated lithium ion 
in oxygen octahedron voids of the anatase ele- 
mentary cell. Occupation of positions 1 and 2 is 
equal in the poor Li-anatase phase and different 
in the rich Li-titanate phase; their orientations 
are antiparallel in the neighbouring octahedra 
(sublattices k = 1, 2). 
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Hubbard projection operator Xfji, where i is the lattice site index (i.e. the elementary cell in- 
dex), k = 1, 2 is the sublattice index (the octahedron index), and the state p = corresponds 
to an empty octahedron while the states p = 1, 2 denote the lithium ion in positions 1 or 2, re- 
spectively. The Hubbard operator formalism reflects the microscopic structure of the system and 
is very convenient for further calculations. 

Alternatively, in the pseudospin formalism localization of lithium in a certain position can be 
described by the pseudospin operator Si k = ( — 1) H^/fe — X 22 ). However, one should also take 
into account the total occupation of the void = X^ + X 22 . Such an approach separates dipole- 
dipole (pseudospin) and particle-particle interactions in the spirit of the Blume-Emery-Grifflths 
(BEG) model [ED]. 

Finally, symmetrized linear combination of the averages 

n± = -(ni±n 2 ), r,± = ^(si±s 2 ) (n k = {n ik ), s k = (s ik )) (1) 

inherit symmetry properties of the system, thus being the order parameters of possible phase 
transitions (see appendix). Namely, n+ (it transforms according to the irreducible representation 
(IR) A\ g of the point group D^h) corresponds to the lithium concentration (the average occupation 
of octahedral voids), n_ (IR B 2u ) is the difference of the void occupations in the sublattices 1 and 
2, ry+ (IR A 2u ) is the polarization along the z axis, and ?y_ (IR B\ g which corresponds to the 
phase transition into the Li-titanate point symmetry subgroup D 2 h) simultaneously describes two 
phenomena: the antipolarization along the z axis (unlike the true antiferroelectric ordering with 
doubling of the unit cell, a mutual compensation of sublattice polarizations occurs here just as 
in the Mitsui model) and the deformation U xx — U yy in the ab plane. Such a coexistence of the 
antipolar ordering and the deformation belonging to the same irreducible representation and, thus, 
described by a common order parameter is called internal piezoeffect. 



3. Phase equilibrium in the framework of the Landau expansion 

The symmetry analysis performed in the previous section can serve as a background for a 
qualitative description of thermodynamics of the considered system in the framework of Landau 
expansion 

F = F + l -ap 2 + hp 3 + hzp A + ^Arj 2 + -Br? - a V - pp , (2) 

where p describes the intercalant (lithium) concentration and hence corresponds to the n + intro- 
duced above, p is the chemical potential, the order parameter rj is proportional to the ry_, and 
the conjugated "field" a describes the applied stress; the expansion coefficients should satisfy the 
following conditions 

A = A + Axp, B > 0, c> 0. 
The equilibrium state of the system is achieved at the minimum of the free energy 

dF 1 

— = ap + bp 2 + cp 3 + -Am 2 -H = Q, (3) 

— = A V + Br, 3 - a = 0. (4) 
or, 

Further calculations are limited to the case of zero "field": a — 0. Then, equation (J4j) can have 
either trivial solution r, = or nonzero one: 



A + Bt, 2 = => 770 = ±y/-(A Q + A lP )/B. (5) 

At B > 0, r,Q takes on a real value under the condition Aq + A\p < which means Aq > and 
A\ < giving 

2 \Ai\p-Ap A 

% = — B ' P> \M[ (6) 

Let us consider the cases r, — and r, ^ separately. 
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I. 77 = 

As follows from equation J3]) 



p(p)=M> (p(p) = ap + bp 2 + cp 3 . (7) 



This equation could have three solutions in a certain region of chemical potential values (i.e. 
the possibility of a phase transition with the jump of p) if the extrema of the function ip(p) 
exist, i.e. the equation 

fitni n\ „ 

= a + 2bp + 3cp 2 = 



dp 

has nonzero solutions 

1 (8) 



1 

Pl,2 = 7T 

3c 



-b ± \Jb 2 — 3ac , 
which imposes a condition on the Landau expansion coefficients 

b 2 - 3ac > 0. (9) 
An equivalent condition could be obtained by setting the second derivative to zero 

-Pf- = 2b + 6cp = 0, 
op'' 

which gives the ordinate of the inflection point 

p* = -&/3c, (10) 

and demanding a negative value of the first derivative at this point d<p(p)/dp\ p=p * < 0. Since 
the curve ip(p) always crosses the inflection point (which is the symmetry centre of the curve), 
this point is also crossed by the line of the phase transition occurring at the following value 
of chemical potential 



3c V 9c 

Considering that the parameter p describes concentration, we have an additional condition 
p ^ 0: both solutions of equation ([8]) are positive if b > 0. 

II. 7] = ±7] 

After the identical calculations one can obtain an expression similar to equation |(7| but 
slightly renormalized: 

, 2 3 \Ai\ 2 A \Ai\ . , 

ap + bp +cp = p,, a = a ~~2B~-> 2B~ ' ^ 

The phase transition exists if 

b 2 - 3ac> (13) 
and it occurs at the following value of chemical potential 

* b fib 2 \ , , 



As follows from the above considerations, the behaviour of the system at change of the chemical 
potential depends on the values of Landau expansion coefficients. So, further considerations are 
limited to the case which qualitatively describes the phase transition between the poor phase I 
(r/ = 0, p — > 0) and the rich phase II (rj ^ 0, p — * 0.5) in the lithiated anatase. 

Comparing conditions |(9j) and (fT3|) one can derive the condition 

3 \A I 2 

3ac- y-c < b 2 < 3ac (15) 
2 B 
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describing the case when only one solution for p exists in the phase I while there are three possible 
solutions in the phase II at nonzero t)q. Combining the equations for equilibrium values of p and 
the condition of the first order phase transition Fj = Fu between the phases I and II, we obtain a 
set of equations for values of p\ , p\\ and p at the phase transition point 



api + bpf + cpf = p, 

apii + bpfj + cpfj = p, 

1 2 1. a 1 4 

-ap l + -bp 1 + -cp 1 - ppi = 

1 2 



(16) 



1— 
0) 

c 
<a 

0) 
CD 
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7 = 
■ - 77*0 
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A qualitative illustration of such a phase 
transition, which corresponds to the case in the 
lithium intercalated anatase, is given in figured] 
for the following set of Landau coefficients (O : 
a = 10.8, b = -35.2, c = 45.7, A = 0.2, 
A\ = —9.1, 5 = 9, This set satisfies the con- 
dition lfT5]) and, as follows from the free energy 
analysis, at the rise of chemical potential the 
first order phase transition from the branch I 
to the branch II (fT6|) precedes the possible sec- 
ond order phase transition with a continuous 
growth of the nonzero value of p at p(p) = p c . 

Due to the dependence on parameter p of 
the coefficients of parameter r\ both of them 
have simultaneous jumps at the point of the 
first order phase transition (figure [3]) . Parame- 
ter T] should be considered as a true order pa- 
rameter of this transition because it is exactly 
equal to zero in the initial phase. It should be 
also noted that (due to proximity of the phase transition points) at increase of chemical poten- 
tial the system could "pass through" the first order phase transition point (a metastable state) 
and the second order phase transition does occur followed by the first order phase transition in 
the extremum point of p{p). All the above considerations are valid for the case p — const. How- 
ever, the lithiated anatase corresponds rather to the system with the fixed lithium concentration 
(p = const). In this case the system separates into phases with concentrations p\ and p\\ (figure [3]) 
and respective weights w\ and u>n, so pfi xe d = W \P\ + ^iiPii- The chemical potential of the system 
is constant and is equal to the chemical potential value at which the first order phase transition 
occurs in the p = const regime. 
Using the equality 

d 2 F , on 2 da 
— = A + Wr 1 2 = — , 
or] z or] 

one can calculate the susceptibility which describes the reaction of the order parameter r\ (the 
deformation U xx — U yy ) with respect to the "field" a (this susceptibility is related to the elastic 
modulus of the system) 



Figure 2. Dependence of free energy on chemical 
potential for branches with zero and nonzero or- 
der parameter r\: the first order phase transition 
takes place. 



X 



dT1 - \A 



lA^p + SBr, 2 



Its explicit form depends on the phase 



r/ = 0: X i = [4)-|^lM \ 

V = Vo- Xii = 2 Pi Iwi - AiP 1 ■ 



(17) 

(18) 
(19) 
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Figure 3. A jump of order parameters p (left) and r\ (right) at the first order phase transition. 
Thermodynamically stable solutions are marked with the solid curve. In the p = const regime 
the phase separation into phases with p\ and pu takes place. 



Due to the phase transition, the susceptibility has a jump whose value can be calculated from 
expressions lfl8l) and (fl9| using solutions of the set (fl6l) . 

Utilizing the equality Aq — \A±\p = —Br/ 2 , one can rewrite the susceptibility lfT7|) at rj 7^ in 
the form 

X=[2BrfY 1 . (20) 

Thus, the susceptibility could diverge in both cases lfT8]) and ([201 approaching the point of the 
possible second order phase transition from the respective direction. 

4. Lattice model for the lithium intercalated anatase 

4.1. Model Hamiltonian and thermodynamics in the mean field approximation 

As demonstrated in the previous section, the Landau expansion combined with the symmetry 
analysis gives a good qualitative picture of the phase separation in the lithiated anatase. However, 
a detailed description of temperature dependent thermodynamic properties of the system could be 
derived only in the framework of a microscopic approach. Let us construct a model Hamiltonian 
of the lattice gas type: 

H = H 1 + fli nt + fldef) (21) 

where besides the single particle term Hi the interparticle interaction H mt and the lattice defor- 
mation fldef are taken into account: 

& = EEE( £ » " ») X ™ - h^sa+sa) - A J>i- Si2 ), 

i k p i i 

i^j kl pq 

H de( =^NCU 2 7 (22) 

where 

A = aU, U = U xx -U yy ; fc,J = l,2, p,g=0,l,2; 

fi is chemical potential of the intercalant particles (let us assign e as its origin) , h is the external 
electric field, A is the deformation induced internal field, U is the effective deformation in the 
XY plane (a = d(Cl/N)/dU defines the stress), W^(i,j) are the interaction energies between the 
particles in the respective positions. Hence, we take into account both the semiphenomenological 
deformational shift of the lattice energy fl^f and the effective internal staggered field A (similar 
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to the one in the Mitsui model) which appears due to the intercalation induced lattice deformation 
making the lithium intercalation positions inequivalent in the pairs. Unlike the ordinary lattice gas 
approach, the proposed model considers two equilibrium positions for intercalated particles. Such 
multistate models (see, e.g. [ [4TJ [42]) are rather rare since monopositional intercalated materials 
are the most common. Deformational effects have been also taken into account because an effective 
potential for the Li ions was affected by the intercalation- induced distortion of the host [ 43j . 
In the mean field approximation, Hamiltonian l(2T|) becomes linear 



H MFA = NE + J2J2J2 H *P X ik- (23) 

i k p 

Taking into account the symmetry properties of the Fourier transforms of the interaction energies 
in the centre of the Brillouin zone 

wll = wll = w\l = wil w\l = wll = wg = wg, 

w\l = wl\ = wll = wll wll = wll = wll = wll, 

one can write down the expressions for average values and effective fields 

e = \wll ((x±y + (xir + (x 1 , 1 ) 2 + (xff) 

+ Wll((X?)(X?) + (X?)(X?)) 
+ Wll{{X?){X?) + {X?){X?)) 
+ WH((X±i)(X?) + (X?)(X?)), 
Hix = -n - h - A 

- (WH{X?) + WH{X?) + Wg{X?) + WH{Xf)) , 
H 12 = -fJt + h + A 

- (wll(x?) + wlKx?) + w£(x?) + wll(x?)) , 

H21 = -fJt + h - A 

- (w$(x?) + WH(X?) + W$(X?) + wll(xf)) , 
H 22 = -fi-h + A 

(Wll{Xf) + Wg(X?) + WU{Xf) + WIKXI 1 )) , (24) 
Average occupations of the positions can be obtained as solutions of the selfconsistency equation 

set 

(Xf> = Ve"^, (25) 
where partition functions of sublattices are as follows: 

Z k = 1 + e-P H " + e~P Hk2 (26) 

and thermodynamically stable solutions are chosen according to the criterion of the minimum of 
grand canonical potential 

Sl/N = E + \CU 2 -ein(ZiZ 2 ). (27) 
For the (n, s)-representation (introduced in section [2]) selfconsistency equations look like 

^ = Z -i(e-^ +e -^), 

s k = {-if^Z^ 1 (e-^ 1 - e^* 2 ) , (28) 
where the term Eq and the effective fields are as follows 

Eo = \ [W+{n\ + n\) + W^{s\ + s\) + 2W+n x n 2 - 2W^ Sl s 2 ] , 
Hu = -ii - h - A - ± (W+ ni + W nSl + W+n 2 - W u s 2 ) , 
H 12 = -(i + h + A-± (W+ ni - W nSl + W+n 2 + W u s 2 ) , 
H 21 = -tx + h-A-\ {W+ ni + W uSl + W+n 2 - W n s 2 ) , 
H 22 = -n - h + A - \ (W&rn ~ W uSl + W+n 2 + W n s 2 ) , (29) 
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and new combinations of interaction energies are introduced 

W± = W\\ ± W\l , W± = W\\ ± W\l . (30) 

Finally, the symmetrized combinations of averages ([T]) (which can be the order parameters of 
the system) are found from the following set 

n± = \ [Z^ 1 (e-^ + e-^) ± Z 2 X (e^ 21 + e^"™)] , 
n± = \ [Zf 1 (e-^ u - e"^ 12 ) =F Z 2 - x (e-^- - e"^)] , (31) 
with the respective definitions for E$ and the effective fields 

E Q = \ [W++n\ + W+-n 2 _ + W—T]\ + W-+r] 2 _] , 
ffn = -M - h - A - \ {W++n + + W+-n^ + W__r? + + W^+rj-) , 
H12 = -^ + /i + A- i {W ++ n+ + W+-n^ - W—ri+ - W- + Tj^) , 
#21 = -n + h - A - \ {W++n + - W+_n_ - 77+ + VK_+r;_) , 
#22 = -M - ft + A - | (iy++n + - W + _n_ + VK__ 77+ - VK_+r;_) , (32) 

as well as symmetrized interaction energies 

W+± = ± W- ± = Wxx ± W^. (33) 

Considering the definition a = d(Vt/N)/dU with account of the equality A = aU and the 
expression for grand canonical potential (|27l) . the equation for the deformation is obtained 

U=^(a + r 1 .), (34) 

where a = a /2a is a scaled dimensionless stress. As it follows 

ICU 2 = k A {a + r,_) 2 , A = fc A (& + r?_), (35) 

where fc A = 2a 2 /C. Thus, the deformation lj34|) can occur spontaneously (giving rise to the anti- 
symmetrical internal field A) due to the appearance of the order parameter r]_ even at the absence 
of the stress. 

The deformation U is a proper variable for the grand canonical potential fl i|27p but in our case 
it is convenient to deal with the stress a (conjugated to U). Performing the Legendre transformation 

dO = • ■ • + adU = ■■■ + d(aU) - Uda, 

one can build the desired form of grand canonical potential O 

n = n-aU. (36) 

The respective deformational term of the thermodynamic potential looks like 

^CU 2 -all = k A { V 2 _~a 2 ). 

4.2. Phase diagram of the ground state 

At zero temperature, the homogenious system, which is described by Hamiltonian (|2"Tj) . could 
reside in one of the nine possible states \pxP2} (let us also use a more descriptive notation where 
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"up" and "down" arrows indicate occupied positions, e.g. 1 10) = |f0)) with the following values of 
thermodynamic potential (|36j) : 



|00) EE 


|0): 


n = 


-k A a 2 , 


|T0) EE 




hx = 


- f i-h-W 1 -k A (a+^) 2 


U0> EEEE 


|2): 


n 2 = 


-fl + h-Wi -fc A (cr- ±) 2 


\0i) = 


|3): 


n 3 = 


-fi + h-w 1 -k A (a + ±) 2 


|0T) = 


|4): 


Qi = 


-li- h- Wi -k A {a- ±) 2 


ITT) = 


|5): 


n 5 = 


-2(i-2h-2W 2 - k A a 2 , 


|U) = 


|6): 


n 6 = 


-2^i + 2h-2W 2 -k A a 2 , 


in) = 


|7>: 


n 7 = 


-2fi-2W 3 -k A {a + l) 2 : 


ITT) = 


|8): 


n 8 = 


-2n-2W 3 -k A (a-l) 2 ; 



(37) 



where 



W 2 
W 3 



= 1(^4 



W-) = \{w±t + wl 2 ), 



(38) 



As is obvious from expressions (|37)) . at a > (in particular at a — > +s) the levels 1, 3, and 7 
always lie below the levels 4, 2 and 8, respectively. So, the latter will not be further considered. 

Setting equal thermodynamic potentials of different phases, one can obtain a set of equations 
for the respective phase transitions: 

-h-Wi-kA(v+\), 
h-Wx-k^a + l), 

-h-W 2 , 
h-W 2 , 

-W 3 -k A (a+i), 

0, 
0, 

(W 3 -W 2 ) + k A (a+±), 
-(W 3 -W 2 )-k A (a+i), 
-h+{W 1 -2W 2 ) + k A (a+\), 
h+(W 1 -2W 2 )+k A (a+\), 
h+{W 1 -2W 3 )k A {a+\\ 
-h+{W 1 -2W 3 )k A {d+\), 
ih+{W 1 -2W 2 ) + k A {a+\), 

-3h+(W 1 -2W 2 ) + k A (a + \). (39) 

Various possible phase diagrams of the ground state in the [i-h plane (which are calculated from 
equations (|39|) ) are depicted in figures 2HZ] using the following designations: 



0) <-> 


1) 




0) <-> 


3) 




0) <-> 


5) 




0) 


6) 




0) 


7) 




1>~ 


3) 


/i 


5) <-> 


6) 


/i 


5) 


7) 




6) 


7) 




1)- 


5) 




3) <-> 


6) 




1>~ 


7) 




3) <-> 


7) 




1)- 


6) 




3) 


5) 





I). 



/i = -W 3 - k A (a 
Hi = -Wi-k A (a+\), 
» 2 = W 1 -2W 3 -k A (a+Z), 
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Figure 4. The most general form of the ground 
state phase diagram: existence of intermedi- 
ate half-filled phases and the possibility of a 
direct transition from the empty state to the 
full one. 



Figure 5. The ground state phase diagram 
with the phase transition between the empty 
and half-filled states which corresponds to 
the Li-poor - Li-rich phase equilibrium in the 
lithiated anatase. 



H3 = Wi-2W 2 + k±{*+\) ) 

h 1 = W 3 -W 2 + k A {a + ±). (40) 

It is evident that the form of the phase diagram depends on the values of model parameters, so we 
shall analyse some important cases below. 

In the most general case (figure 0]), the phase diagram consists of empty, half-filled, and full 
states which differ in occupation and polarization of sublattices. Two-sublattice nature of the model 
demonstrates itself in the intermediate half-filled states and the central nonpolar full state thus 
being noticeably different from the respective ground state diagram of the BEG model. It should 
be stressed that due to the exclusion of the states 2, 4, and 8 the phase |t/0) = |1) is symmetrical 
to the phase |0|) = |3). In the considered case, both transitions between the empty and half-filled 
phases as well as between the empty and full phases are possible. The latter transition takes place, 
e.g. at zero field h which is contrary to the situation in the lithiated anatase where only half-filled 
phases are accessible. 

Coexistence of the Li-poor and Li-rich phases in the lithiated anatase is fairly described by 
a phase diagram in figure [H For any value of the external electric field h (in particular, at zero 
field) the system can pass from the empty state to the half-filled state only. In the regime of fixed 
concentration such a phase transition manifests itself as a phase separation. 

Such a phase diagram exists if conditions [i\ < [i 2 < are satisfied. Taking into account 
definitions (|4*0|) and f38|) and neglecting deformational effects one can rewrite the above inequality 
in a more clear form 

W\l < W\\ < 0. (41) 

Thus, this type of a phase diagram exists if the interactions between the orientational states in 
different sublattices are repulsive with the interaction between the unlike states (with different p, 
i.e. I TO) and |0f) or |J,0) and |0|}) being stronger than the interaction between the like states (with 
identical p, i.e. ||0) and |0J,) or ||0) and |0T))- These conclusions are in full qualitative agreement 
with the predictions of the molecular dynamics modelling [ [37] . 

In view of deformational effects, the inequalities ^1 < ^2 < transform to 

- k A (a + |) < -WU - k A (& + |) < -Wg + k A (a + \). (42) 

Thus, application of the stress a favours the phase diagram in figure[5]because /Zi and ^2 are shifted 
to the left ~ a and ^ is moved towards the right ~ a. So the domains of the phases ||0) and |0J.) 
expand. However, at h — the difference \i 2 — fJ-i = — \k A does not depend on the stress 
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Figure 6. The effect of stress imposition: the ground state phase diagram in the /i-ct coordinates 
at h = (left) and a family of ground state phase diagrams for different values of the stress 
(right). Other model parameters have the following values: W++ = 1, W-\ — = 1.5, W-- = 0.2, 
W-+ = 0.3, k A = 0.05. 
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Figure 7. A reduced form of the ground state 
phase diagram: only the empty and full phases 
coexist. 



(as one can see on the phase diagram [i-a in 
figure [6]) and half-filled phases exist if W{\ < 
— 2&a- A strong enough stress a can "open" 
these phases even if conditions l|4Tj) are not 
satisfied. The above rationales are confirmed 
by the family of phase diagrams in figure [6] 
calculated for different stress values <r (here- 
inafter all model parameters are given in the 
dimensionless units normalized by W ++ = 1; 
in the considered case =0.75, Wjf = 0.5, 
W\\ — —0.1, Wjf = —0.15 in such a manner 
satisfying conditions (|4"Tj) ). 

If conditions (J32J) fail, the half-filled phases 
are suppressed (figure[7J). In addition, if hi ^ 0, 
then the central full nonpolar phase vanishes 
and the ground state phase diagram reduces to 
the respective one for the BEG model. 



4.3. Temperature behaviour of phase diagrams and phase separation 

Since the considered model inherits features of the Mitsui and BEG models, one can expect a 
quite complex thermodynamical behaviour and the above analysis of the ground state has proved 
these anticipations. As it is seen from the temperature axis complemented analogue (figure[8]) of the 
diagram in figure [H the lines of phase transitions form the surfaces and some new phases appear. 
A comprehensive analysis of the obtained diagram is too complicated and goes beyond the scope 
of this research, so we consider the case with intermediate half-filled phases corresponding to the 
intercalated anatase. 

With this in mind let us analyse the phase diagram at the absence of the external electric field 
h (figure [9]). At low temperature, in full agreement with the ground state diagram, there are three 
phases: "empty", "half-filled", and "full" (due to the temperature "blurring" at higher temperatures 
these names loose their exact meaning). At high temperatures in the whole range of chemical 
potential there is only "empty" phase where all order parameters are zero except the concentration 
n + which changes monotonously. 
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Figure 8. The three-dimensional phase diagram in the fi,~h~@ coordinates. Model parameters 
have the following values: W++ = 1, W+- = 1.5, W— = 0.2, W-+ = 0.3, k A = 0.05, a = 0. 

The "half-filled" phase provides the most 
complex behaviour of the order parameters. As 
is seen from the "cross-section" at the chemi- 
cal potential value /z = —0.38 (figure fTOj) . the 
rise of temperature leads to the suppression 
of the dipole-dipole ordering (order parame- 
ters 77+ and 7?_) and it completely vanishes at 
the line of the second order phase transition 
which is located inside the "half-filled" phase. 
The phase itself is limited from the top side by 
the line of the phase transition with zeroing of 
n_ . The upper part of the phase is separated 
by the tricritical points (TCPs) marking the 
change of the phase transition order from the 
second to the first one. The "full" phase is also 
limited from the top side by the line of the sec- 
ond order transitions where r?_ — > 0. 

Behaviour of order parameters at the 
change of chemical potential (figure fTOj) clearly 
distinguishes the phases separated by the lines 
of the first order phase transitions (e.g. a char- 
acteristic feature of the "half-filled" phase is 
n_ ^ 0). In a certain temperature range inside 
the "half-filled" phase, the above mentioned 
phase transition between the polar and non- 
polar states takes place. 
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Figure 9. The phase diagram /i-O at the absence 
of the electric field. Other model parameters are 
as follows: W++ = 1, W+- = 1.5, W— = 0.2, 
W — |- = 0.3, k/\ — 0.05, <t = 0. Phase names 
("empty", "half-filled", and "full") correspond to 
the ground state. There are two tricritical points 
(TCPs) on the phase transition line limiting the 
"half-filled" phase from the top. 
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Figure 10. Phase diagram "cross-sections": temperature dependences of order parameters in the 
"half-filled" phase (left: /i = —0.38; in the picture scale curves ?7+ and r\- overlap) and depen- 
dences of order parameters on chemical potential exhibiting phase transitions of the second and 
the first orders (right: = 0.136). Other model parameters have the following values: W++ = 1, 
W+- = 1.5, W— = 0.2, W-+ = 0.3, k A = 0.05, a = 0. 



In a wide temperature range, the appearance of the phase diagram "chemical potential fi - 
stress a" (figure ITTj) closely resemble the ground state one (figure [6]). However, the further rise of 
temperature leads to a fast suppression of the "half-filled" phase. 

Since in the regime of a fixed chemical potential {fi — const) the "empty" and "half-filled" 
phases on the phase diagram /i-6 (figure [9]) are separated mainly by the line of the first order 
phase transitions, the system separates into "poor" and "rich" phases (figure fl2|) in the regime of 
the fixed concentration (n + = const). As one can see, in a wide region of low temperatures a 
separation into concentrations n + = and n + = 0.5 occurs which well reproduces the coexistence 
of Li-poor and Li-rich phases in the intercalated anatase. The phase separation region narrows at 
heating and finally closes up at the point corresponding to the tricritical point in figure [9l Another 
short line of phase transitions (starting at the kink of the right boundary) relates to the suppression 
of polar states in a "half-filled" phase. 




Figure 11. A family of phase diagrams "chemi- 
cal potential fi - stress a" for various tempera- 
ture values at h = 0. Other model parameters 
have the following values: W++ = 1, W-\ — = 
1.5, = 0.2, W-+ = 0.3, k A = 0.05. 



Figure 12. The diagram of the phase separa- 
tion into "poor" (n+ = 0) and "rich" (n+ = 

0. 5) phases in the regime n+ = const (W++ = 

1, W+_ = 1.5, W— = 0.2, W_+ = 0.3, 
k A = 0.05, b = 0). 
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A question arises why in the microscopic 
model we should deal with four order param- 
eters while the Landau expansion is quite suc- 
cessful with only two? First of all, the semiphe- 
nomenological description is qualitative only 
and it just demonstrates the possibility of the 
first order phase transition with simultane- 
ous jumps of concentration and the order pa- 
rameter r\ as the minimal set necessary to 
describe the phase separation in the litiated 
anatase. Such a picture corresponds to a direct 
phase transition between the "empty" and "full" 
phases (see the phase diagram at zero external 
field in figure [9]) when other order parameters 
(n_ and r] + ) are always equal to zero. The de- 
pendences of the "active" order parameters n + 
and r\- on chemical potential (figure fl4|) closely 
resemble the respective curves for p and T] ob- 
tained by the Landau expansion (figure [3]). But 

if one should take into account the "half-filled" phase (what is inevitable for description of the lithi- 
ated anatase), the order parameters (n_ and r/ + ) became nonzero and the full set of four parameters 
should be considered as it has been done above. 
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Figure 13. The phase diagram fj,-Q at h = 
for the case of the first order phase transition 
between "empty" and "full" phases (n_ = and 
r] + = in both of them; W++ = 1, W+- = 0, 
W— = 0, W-+ = 0.3, k A = 0.05, a = 0). 
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Figure 14. A jump of order parameters n+ (left) and 7/_ (right) at the first order phase transition 
illustrating the previous phase diagram. Thermodynamically stable solutions are marked with 
the solid curve. Other model parameters have the following values: W++ = 1, W-\ — = 0, W-- = 
0, W-+ = 0.3, k A = 0.05, = 0.13, h = 0, o = 0. 



5. Conclusions 

To sum up, the present study was inspired by two features of the lithium intercalated anatase: 
coexistence of poor and rich phases and two possible localizations of Li ion in the oxygen octahedron 
along the c axis. The second one implies the possibility of (anti)polar ordering what is beyond the 
scope of a simple lattice gas model well describing a phase separation. So a model of Blume-Emery- 
Griffits type has been used which takes into account both particle-particle and dipole-dipole terms. 

Another less obvious peculiarity of the considered compound is the intercalation induced de- 
formation of lattice: Li-rich phase has a lower symmetry (the axes a and b become inequivalent) 
resulting in preferential occupation of one position of the mentioned pairs which has an alternating 
orientation in the neighbour octahedra (i.e. in different sublattices). Performed symmetry analysis 
explains this phenomena by the possibility of internal piezoeffect: the deformation in the ab plane 
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as well as the appearance of an effective internal staggered field (causing the ordering of antiferro- 
electric type like as in the Mitsui model) both belong to the same irreducible representation of the 
initial high-symmetry anatase phase and, hence, are described by the common order parameter. 
Thus, increase of the intercalant content could result in a phase transition with simultaneous jumps 
of the average occupation and antipolarization (the latter accompanied by the jump of dielectric 
susceptibility) as it has been proved by the Landau expansion. 

The microscopic approach, combining the abovementioned features of both the BEG and Mitsui 
models, gives semiquantitative description of phase coexistence in the lithiated anatase. Analysis 
of the ground state phase diagram confirms a possibility of the phase transition between "empty" 
and "half-filled" phases which corresponds to the phase separation into Li-poor and Li-rich phases 
in the crystal. As the model predicts, such a separation remains near constant in a wide tem- 
perature range. The microscopic approach could easily reproduce the Landau expansion results 
as the particular case of the "empty"-"full" transition described by the two order parameters. But 
the presence of the "half-filled" phase makes it necessary to deal with the full set of the order 
parameters allowed by the crystal symmetry. 

However, some issues are still open. The model predicts that external stress should shift the 
phase transition between empty and half-filled phases to the lower values of chemical potential. An 
experimental evidence of this conjecture is still missing. The real average occupation in the Li-rich 
phase is 0.55-0.6 instead of the value n + = 0.5 in the half-filled model phase. This deviation could be 
explained by a multidomain nature of the Li-rich phase containing "impurities" of the full-occupied 
LiTiC>2 phase while the model phases are monodomain by definition. The same explanation applies 
to the issue of absence of the total polarization in the half-filled phases: at zero external field 
these phases with opposite polarizations could coexist in different domains providing a full mutual 
compensation. A similar mechanism of the mutual compensation of polarization in Li-"chains" with 
an opposite Li orientation is supported by the molecular dynamics simulations [ 37J. 
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Appendix: Table of the point group D± h symmetry transformations for the 
symmetrized averages 



Table 1. Transformations of the symmetrized linear combinations of the averages which corre- 
spond to the irreducible representations (IR) of the point symmetry group D^. The operations, 
which also belong to the lower symmetry subgroup D2h, are marked by asterisk; they keep the 
symmetrized combination rj- (IR B\ g ) invariant. 
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